Nonlinear guided waves and spatial solitons in a periodic layered medium 
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We overview the properties of nonlinear guided waves and (bright and dark) spatial optical soli- 
tons in a periodic medium created by a sequence of linear and nonlinear layers. First, we consider a 
single layer with a cubic nonlinear response (a nonlinear waveguide) embedded into a periodic lay- 
ered linear medium, and describe nonlinear localized modes (guided waves and Bragg-like localized 
gap modes) and their stability. Then, we study modulational instability as well as the existence 
and stability of discrete spatial solitons in a periodic array of identical nonlinear layers, a one- 
dimensional nonlinear photonic crystal. Both similarities and differences with the models described 
by the discrete nonlinear Schrodinger equation (derived in the tight-binding approximation) and 
coupled-mode theory (valid for the shallow periodic modulations) are emphasized. 



I. INTRODUCTION 

Spatially localized waves (or intrinsic localized modes) 
in nonlinear lattices have been an active research topic 
during last years. In application to the problems of non- 
linear optics, such modes are known as discrete spatial 
solitons, and they have been described theoretically [QJ|] 
and recently observed experimentally |^] in periodic ar- 
rays of nonlinear single-mode optical waveguides. A stan- 
dard approach in the study of the discrete spatial solitons 
in optical super-lattices is to employ the properties of 
an effective discrete nonlinear Schrodinger (NLS) equa- 
tion, which can be derived under some assumptions, sim- 
ilar to the tight-binding approximation in the solid-state 
physics [Q. On the other hand, nonlinear localized waves 
in a system with a weakly modulated optical refractive 
index are known as gap solitons [|J . The similar problems 
and methods of their solutions appear in other fields, such 
as the nonlinear dynamics of the Bose-Einstein conden- 
sates in optical lattices || . 

However, real experiments in guided-wave optics are 
conducted in the structures of more complicated geome- 
tries and, therefore, the applicability of the tight-binding 
approximation and the corresponding discrete equations 
become questionable. Moreover, one of the main features 
of the wave propagation in periodic structures (which 
follows from the Floquet-Bloch theory) is the existence 
of forbidden transmission band gaps and, therefore, the 
nonlinearly-induced wave localization can be also possi- 
ble in the form of the so-called gap solitons located in each 
of these gaps. However, the effective discrete equations 
derived in the tight-binding approximation describe only 
one transmission band surrounded by two semi-infinite 
band gaps and, therefore, the real fine structure of the 
band-gap spectrum associated with the wave transmis- 
sion in a periodic medium is lost. On the other hand, 
the coupled-mode theory of the gap solitons Q describes 
only the modes localized in an isolated narrow gap, and it 
does not allow to consider simultaneously the gap modes 
and conventional guided waves localized due to the total 
internal reflection. However, the complete band structure 
of the transmission spectrum and simultaneous existence 



of localized modes of different types is very important in 
the analysis of the stability of nonlinear localized modes. 
Such an analysis is especially important for the theory of 
nonlinear localized modes and nonlinear waveguides in 
realistic models of nonlinear photonic crystals (see, e.g., 
the recent paper ||, and references therein). 

In this paper, we consider a simple model of a non- 
linear periodic layered medium where the optical super- 
lattice is formed by a periodic sequence of two linear lay- 
ers of different dielectric susceptibilities, and nonlinear 
waveguides are described by thin-film layers embedded 
into it (?J (see also Ref. ||). In such a case, the effects of 
the linear periodicity and band-gap spectrum structure 
are taken into account explicitly, and nonlinearity enters 
the corresponding matching conditions only, allowing a 
direct analytical study. 

Here, we consider two problems of this kind. In the 
first case (Sec. ||), we study the nonlinear guided waves 
supported by a thin isolated nonlinear layer with the 
Kerr-type nonlinear response, which is embedded into 
a linear periodic structure composed of two layers with 
different linear dielectric constants. We consider the case 
when the layer describes an optical waveguide that sup- 
ports guided waves in a homogeneous medium of the av- 
eraged susceptibility, and assume that the cubic nonlin- 
ear response may be either positive (self-focusing) or neg- 
ative (self-defocusing) . We describe two different types 
of nonlinear localized modes supported by this structure, 
including the analysis of the mode stability. In the second 
case (Sec. Ill), we analyze the nonlinear localized modes 
in an infinite structure consisting of a periodic sequence 
of nonlinear layers, similar to the geometry of the exper- 
iments with discrete optical solitons Q. In this case, we 
study first modulational instability in both self-focusing 
and self-defocusing regimes, and then discuss the prop- 
erties of different types of nonlinear localized modes such 
as bright, dark and "twisted" spatial solitons. We em- 
phasize both similarities and differences with the models 
described by the discrete NLS equation and the contin- 
uum coupled-mode theory. 
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II. NONLINEAR GUIDED MODES 
IN A PERIODIC MEDIUM 



where (3 is the propagation constant, and the real func- 
tion u{x; (3) satisfies the stationary nonlinear equation: 



A. Model 

We consider the wave propagation along the Z- 
direction of a slab waveguide consisting of a sequence of 
linear and nonlinear layers. The evolution for the com- 
plex electric field envelope E(X, Z) in the waveguide is 
governed by the NLS equation, 



dE 
'dZ 



D 



d 2 E 
dX 2 



■e{X)E + g{X)\E\ 2 E = 0, 



(1) 



where D is the diffraction coefficient (D > 0) . The phase 
velocity of the fundamental mode is defined by the func- 
tion s(X), and g(X) characterizes the Kerr-type nonlin- 
ear response. We assume that either the function e{X) 
or g(X) (or both of them) is periodic in X, i.e. it de- 
scribes the periodic layered structure similar to the so- 
called transverse Bragg waveguides created by the nonlin- 
ear thin-film multilayer structures P,|l0|l or the impurity 
band in a deep photonic band gap [[lip 

In order to reduce the number of physical param- 
eters, we normalize Eqs. (|l]) as follows: E(X, Z) = 
ip(x, z)Eqc 1 eZ , where e is the mean value of the function 
s(X), x = X/d and z — ZD/d 2 are the dimensionless co- 
ordinates, d and Eq are the characteristic transverse scale 
and field amplitude, respectively. Then, the normalized 
nonlinear equation has the form 



.dip <9 2 ?/> 



dz 



(2) 



where the real function T{I;x) = d 2 D~ 1 [e(X) — e + 
g(X) I\E Q \ 2 ) describes both nonlinear and periodic prop- 
erties of the layered medium, and / = \ip\ 2 is the normal- 
ized local wave intensity. We note that the system (|^) is 
Hamiltonian, and for the spatially localized solutions the 
power, 



P = 



\ip(x, z)\ 2 dx, 



is conserved. 

At this point, it is important to mention that Eq. (j^) 
describes the beam evolution in the framework of the so- 
called parabolic approximation, valid for the waves prop- 
agating mainly along the z direction (see also Ref. Q). 
In other words, the characteristic length of the beam dis- 
tortion due to both diffraction and refraction along the z 
axis should be much larger than the beam width in the 
transverse direction x. This leads to the condition of a 
weakly modulated periodicity, \e(X) — e\ <C |e|. 

We look for stationary localized solutions of the nor- 
malized equation @ in the form 



4>(x,t) =u{x;(3)e ll3z , 



(3u + —r +F{I;x)u = 0. 
dx 



(4) 



We assume that the basic optical super-lattice is linear, 
and nonlinearity appears only through the properties of 
a thin-film waveguide (or an array of such waveguides) 
Then, if the corresponding width of the wave envelope 
is much larger than that of the waveguide, the thin-film 
waveguide can be modeled by a delta-function and, in 
the simplest case of a single nonlinear layer, we can write 
J- (I; x) — v(x)+S(x)G(I), where the function G(I) char- 
acterizes the properties of the nonlinear thin-film waveg- 
uide, and v{x) = v(x + h) describes an effective potential 
of the superlattice with the spatial period h. 

For such a nonlinearity, localized waves can be con- 
structed with the help of certain matching conditions, by 
using the solutions of Eq. (jj) with T{I;x) — v(x), pre- 
sented in the form of the Bloch-type wave functions [^3) . 
Additionally, this approximation can be easily extended 
to the case of a pe riodic array of thin-film nonlinear 
waveguides (see Sec. Ill A| below). 



B. Band-gap spectrum and localized waves 

If the effective periodic potential v(x) is approximated 
by a piecewise-constant function (such a model is known 
as the Kronig-Penney model), the solution can be decom- 
posed into a pair of counter-propagating waves with the 
amplitudes a{x;(3) and b(x;(3), 

Ub (x; 13) = a(x; /3)e^ (x;,3)x + b(x; f3)e+^ x ' f3)x , (5) 



where /i(x; f3) = yjS — v(x) is the local wavenumber. As 
follows from the Floquet-Bloch theory, for a Bloch-wave 
solution the reflection coefficient r(x; (3) = b(x; (3) /a(x; (3) 
is a periodic function, i.e. r(x; (3) = r(x + h;(3), and it 
satisfies the eigenvalue problem: 



T{x;(3) 



1 



ix;(3) 



t(I3) 



1 

r(x;f3) 



(6) 



(3) 



where T(x; (3) is a transfer matrix that describes a change 
of the wave amplitudes {a, 6} after one period (x,x + 
h) It was found that det T = 1 and, therefore, 

two linearly independent solutions of Eq. (^|) correspond 
to a pair of the eigenvalues r and r _1 . Relation r(/3) 
determines a band-gap structure of the superlattice spec- 
trum: the waves are propagating, if |t| = 1, and they 
are localized, if |t| =/= 1. In the latter case, a nonlin- 
ear waveguide can support nonlinear localized waves as 
bound states of Bloch-wave solutions with the asymp- 
totics \u{x — ► ±oo)| — * 0. The wave amplitude at the 
thin-film nonlinear waveguide is determined from the 
continuity condition at x = 0, i.e. I = Py0 + ) = I{0~) 
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and [du/dx}°* =0 _ + G(I )u x=0 = 0. Then, we use Eq. (||) 
to express the latter condition through the superlattice 
characteristics, 



C. Self- focusing nonlinearity 



Go 

where £ = (£ + + £~ 



G(I ) = CG9), 



(7) 



x=0> 



C ± =M ± 



(1 



(l + r±)' 

and "+" and "— " stand for the characteristics on the 
right and left sides of the waveguide, respectively, i.e. 
v(x) = v + (\x\), for x > 0, and v{x) — v~(\x\), for x < 0. 

Relation (R) allows us to identify different types of 
nonlinear localized states. We notice that such localized 
states exist in the so-called waveguiding regime (Go > 0), 
when (((3) > 0. Additionally, spatial localization can oc- 
cur in the anti-waveguiding regime (Go < 0), provided 
C(/3)<0. 

To study linear stability of the localized solutions, we 
consider the evolution of small-amplitude perturbations 
of the localized state presenting the solution in the form 



ip(x,t) = {u{x)+v(x)e lTz +w*(x)e- iY " z y 



,i/3z 



(8) 



and obtain the linear eigenvalue problem for small v(x) 
and w(x), 



■(P + T)v 



dx 2 



■f v{x)v 

+6{x) [G lV + (Gi 



G )w] = 0, 



(9) 



„„ aw 



v(x)w 

+S(x) [Giw + (Gi - G )v] = 0, 

where Gi = Go + IqG'{Iq). The intensity Iq is calculated 
for an unperturbed solution, and the prime stands for the 
spatial derivative. 

We find that the localized eigenmode solutions of 
Eq. (^) exist only for the particular eigenvalues satisfying 
the solvability condition, Y(-/) = |15[, where 

Y( 7 ) = [G l - cos + r)] [d - cos - r)] - {G l - G f . 

In general, the eigenmode solutions fall into one of the 
following categories: (i) internal modes with real eigen- 
values describe periodic oscillations ( "breathing" ) of the 
localized state, (ii) instability modes correspond to purely 
imaginary eigenvalues, and (iii) oscillatory instabilities 
can occur when the eigenvalues are complex. 

To demonstrate the basic stability results, we consider 
a localized waveguide possessing a cubic nonlinear re- 
sponse, G(I) = a + 7/. Under proper scaling, the abso- 
lute value of the nonlinear coefficient 7 can be normalized 
to unity, so that 7 = +1 corresponds to self-focusing and 
7 = — 1 to self-defocusing nonlinearity. We assume that 
at small intensities the waveguide has a higher refractive 
index, i.e. a > 0. Therefore, below we consider two 
qualitatively different examples corresponding to differ- 
ent signs of the nonlinear coefficient 7. 



We consider the properties of spatially localized waves 
supported by a thin-film waveguide with a self-focusing 
nonlinearity (7 = +1). Such waves do exist in the lin- 
ear limit, and the corresponding propagation constant 
f3(, is defined by the equation C(/?fc) = a - Such localized 
waves correspond to the waveguiding regime (white re- 
gions in Fig. 0, top). We notice that the first band-gap 
(semi-infinite white region at (3 > 19) corresponds to the 
conditions of internal reflection (IR) and, therefore, its 
diffraction properties should be similar to those of the 
conventional guided waves. Localized waves in this band 
resemble the guided waves modulated by a periodic struc- 
ture [see Fig. |l|(b)]. Since an IR wave profile does not 
contain zeros, it is a fundamental eigenmode. Therefore, 
the conditions of the Vakhitov-Kolokolov (VK) stability 
theorem |l6| are satisfied, and the IR states are unstable if 
and only if dP/d(3 < 0. At the critical point, dP/d[3 = 0, 
the linear eigenvalue passes through zero and becomes 
imaginary, as illustrated in Fig. [j] (middle, [3 ~ 27.5). 
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FIG. 1. Top: power vs. propagation constant for the 
nonlinear localized states: solid — stable, dashed — unstable, 
and dotted — oscillatory unstable. Middle: real (dotted) and 
imaginary (solid) parts of the eigenvalues associated with the 
wave instability. Shading marks "waveguiding" (white) and 
"anti-waveguiding" (dotted) localization regimes inside the 
band gaps. Bottom: the localized states corresponding to the 
marked points (a,b) in the top plot; shading marks the areas 
with smaller v. The lattice parameters are h = 1, v(x) = 
for n~l/2<x/h<n, and v (x) = 30 for n < x/h < n+ 1/2, 
where n is integer, and a = 0.5, 7=1. 
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In contrast, band-gaps appear at smaller due to the 
resonant Bragg reflection (BR) by the periodic structure, 
so that the localized waves are similar to gap solitons 
composed of the mutually coupled backward and forward 
propagating waves [see Fig. 0(a)]. For the BR states, 
the VK criterion provides only a necessary condition for 
stability, since the higher-order localized states can also 
exhibit oscillatory instabilities. Indeed, we notice that in 
the linear limit there always exists an internal mode cor- 
responding to a resonant coupling between the BR and 
IR band-gaps, since Y(0^ R ^ — /3p R ^) = 0. We perform 
extensive numerical calculations and find that this mode 
leads to an oscillatory instability of BR waves when the 
value ((3 — Re T) moves outside the band gap; it occurs 
when the wave intensity exceeds a threshold value (see 
Figs. |l|and |). 



exhibit oscillatory instabilities. 




FIG. 2. Example of a resonance that occurs between an 
internal mode of the localized state and a bang-gap edge, 
leading to an oscillatory instability. 



D. Self-defocusing nonlinearity 

Additionally, we study the case of a thin-film linear 
waveguide (a > 0) possessing a self-defocusing nonlin- 
earity (7 = — 1). At small intensities, the waveguide 
supports linear guided waves in the waveguiding regime, 
in both IR and BR gaps, as shown in Fig. |^ (top, white 
regions). Performing the stability analysis, we find that 
IR waves are always stable, while oscillatory instabilities 
appear for higher-order band-gap states. Since at higher 
intensities, i.e. for Iq > \a/j\, the waveguide response 
changes its sign, a new type of localized waves can exist 
in the anti-waveguiding regime (i.e. for the negative ef- 
fective diffraction), bifurcating from the band-gap edge, 
as shown in Fig. || (top, dotted region). It is possible to 
demonstrate that for such waves the VK stability crite- 
rion becomes inverted, i.e. the localized waves are un- 
stable if dP/d/3 > 0. This happens because the signs 
of both the nonlinear response and effective diffraction 
are altered compared to the IR waves supported by a 
self-focusing waveguide. All higher-order band-gap states 



50 



U 

£ 25 
O 
Ph 





1 


i 






















•• a 












FIG. 3. Top: power vs. propagation constant, and (a-c) 
the localized mode profiles for a — 5, 7 = — 1. Notations are 
the same as in Fig. jj. 



III. AN ARRAY OF NONLINEAR WAVEGUIDES 

A. Model and Discrete Equations 

Now we analyze both localized and extended nonlinear 
waves, supported by a periodic array of thin-film waveg- 
uides with the Kerr-type nonlinear response. The cor- 
responding model is a simplified version of the nonlinear 
layered medium usually studied in experiment J3| . In our 
case, we assume again that the nonlinear layers are thin, 
so that the problem can be studied analytically, since 
the nonlinearity enters the periodic matching conditions 
only. 

To simplify our analysis further, we assume that the 
linear periodicity is associated only with the presence of 
the thin-film waveguides, and define the response func- 
tion in the model Eq. (0) as follows 



J- (I; x) = /J (a + r )I)8{x — hn), 



(10) 



where h is the distance between the neighboring waveg- 
uides, and n is integer. As above, the response of thin- 
film waveguides is approximated by the delta-functions, 
and the real parameters a and 7 describe both linear and 
nonlinear properties, respectively. The nonlinear coeffi- 
cient 7 can be normalized to unity, so that 7 = +1 cor- 
responds to self-focusing and 7 = — 1 to self-defocusing 
nonlinearity. The linear coefficient (a > 0) defines the 
linear response (i.e. the response at the lower intensities), 
and it characterizes the corresponding coupling strength 
between the waveguides. 

As has been demonstrated in Ref. Jffl], the stationary 
wave profiles defined by Eqs. (||) and (Q) can be expressed 
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in terms of the wave amplitudes at the nonlinear waveg- 
uides, u n = u(hn), which satisfy a stationary form of the 
discrete NLS equation: 

rjU n + {U n -! + Dn+i) + x\U n \ 2 U n = 0. (11) 

Here U n = \/\£,j\u n are the normalized amplitudes, 
X = sign(C7), and 

i) = —2 cosh(/i/i) + a£, 

n (12) 
£ = smh(V)//i, M = V/3. 



BR IR 




FIG. 4. Characteristic dependencies of the parameters ?7 
and £ on the propagation constant f3. White shading marks 
band gaps. The lattice parameters are h = 0.5 and a — 10. 

Localized solutions of Eq. with exponentially de- 
caying asymptotics can exist for \r)\ > 2, and this con- 
dition defines the spectrum band-gap structure. Char- 
acteristic dependencies of rj and £ vs. the propagation 
constant (3 are presented in Fig. |^, where the band gaps 
are shown as blank stripes. The first (semi-infinite) band 
gap corresponds to the total internal reflection (IR). On 
the other hand, at smaller (3 the spectrum band gaps ap- 
pear due to the resonant Bragg-type reflection (BR) from 
the periodic structure. 



These modes exist in the first transmission band which 
is shifted due to a nonlinear correction to the response 
function. 

One of the main problems associated with the nonlin- 
ear BW modes is their instability to periodic modulations 
of a certain wavelength, known as modulational instabil- 
ity. In order to describe the stability properties of the 
periodic BW solutions, we analyse the linear evolution of 
weak perturbations in the form (||). Now, the functions 
v{x) and w{x) describe the small-amplitude perturbation 
corresponding to the eigenvalue V. Substituting Eq. (||) 
into the original model (g) and ( jlpj ) , we obtain an eigen- 
value problem which has the periodic solutions in the 
form of the Bloch functions, v(x + h) = v(x)e l< * q+KS> and 
w(x + h) = w(x)e l ( q ~ K \ provided the following solvabil- 
ity condition is satisfied, 

[r?(/3 + r) + 2 7 £(/3 + r) Jo + 2 cos(g + K)} 

[ti(0 - T) + 2 7 £(/3 - r)7 + 2 cos(q - K)\ (14) 

= 7 2 £(/3 + r)£(/3-r)/ 2 . 

The spectrum of the possible eigenvalues T is determined 
from the condition that the spatial modulation frequen- 
cies q, which are found from Eq. (|TJ), are all real. 




FIG. 5. Development of modulational instability in a 
self-focusing medium for a slightly perturbed unstaggered BW 
solution with 7o — 0.44. The lattice parameters are a = 3 and 
h = 0.5, and 7 = +1. 



B. Modulational instability 

First, we analyze the properties of the simplest (peri- 
odic) solutions of the model (^) and ( |l0| ) with the equal 
intensities at the nonlinear layers, Iq = \u n \ 2 = const. 
Such Bloch-wave (BW) extended modes are character- 
ized by the wave number in the first Brillouin zone, 
—it < K < 7r, so that u n = uoe lKn and 

7o = -[2cos(if)+r ? ]/( 7 £). (13) 



In what follows, we consider two characteristic cases 
of the stationary BW profiles, when they are (i) un- 
staggered (K = 0) or (ii) staggered (K — tt). Mod- 
ulational instability of the periodic BW solutions has 
been earlier studied for the Bose-Einstein condensates 
in optical lattices |l7| in the mean-field approximation 
using the Gross-Pitaevskii equation, which is mathemat- 
ically equivalent to Eq. (^) with J- (I; x) — v{x) + xl- It 
was demonstrated that the unstaggered solutions are al- 
ways modulationally unstable in a self-focusing medium 
(x = +1), and they are stable in a self-defocusing 
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medium (\ = — 1). It can be proven that the similar 
results are also valid for our model. We note that at 
small intensities the modulational instability in the self- 
focusing case corresponds to long-wave excitations, as il- 
lustrated in Fig. 0. 

The staggered BW modes in a self-defocusing medium 
(x = —1) are always modulationally unstable [jl7). On 
the other hand, the properties of such modes in a self- 
focusing medium (x = +1) are less trivial. We find that, 
in this case, the oscillatory instabilities (i.e. those with 
the complex T) can appear due to resonances between 
the modes of different bands. Such instabilities appear 
in a certain region of the wave intensities in the case of 
shallow modulations, when a is below a certain threshold 
value, as shown in Fig. ^. 




FIG. 6. Modulationally unstable staggered BW modes 
(gray shading) in a self-focusing medium, shown as the inten- 
sity (7o) vs. the lattice parameter a (at h = 0.5, 7 = + 
Dashed line is the analytical approximation ( Jl 5| ) for the 
low-intensity instability threshold. 



which predicts the stability of the staggered modes in a 
self- focusing medium fll| . Numerical and analytical re- 
sults confirm that our solutions are indeed stable in the 
corresponding parameter regions. 
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FIG. 7. Unstable modulation frequencies (gray shading) 
vs. the intensity Io of the staggered BW modes (a = 3, 
h — 0.5). Solid line shows the instability with the largest 
growth rate. 

We find that in a self-focusing medium the staggered 
waves are always stable with respect to low-frequency 
modulations, see Fig. [7]. However, at larger intensities 
unstable frequencies are shifted towards the edge of the 
Brillouin zone, q — it, which has the largest value of the 
growth rate. The corresponding modulational instability 
manifests itself through the development of the period- 
doubling modulations, as shown in Fig. |8|. 



It is interesting to compare these results with the 
results obtained in the framework of the continuum 
coupled-mode theory, valid for the case of a narrow band 
gap, i.e. for small a and small Iq. As was demonstrated 
earlier [Q, in such a case the oscillatory instability ap- 
pears above a certain critical intensity, which is propor- 
tional to the band-gap width. In our case, we find the 
following asymptotic expression for the low-intensity in- 
stability threshold, 

7 4 cr) ~a + 2V2ha 3/2 /TT + 0{a 2 ). (15) 

Since the band-gap width is 2a /h + 0(a 3 / 2 ), we observe 
a perfect agreement with the results of the coupled-mode 
theory (see Fig. [| dashed line). 

In the limit of large a (or large io), the BW solution 
is composed of weakly interacting waves, each localized 
at the individual waveguide. Then, the collective dy- 
namics of the nonlinear modes in such a system can be 
studied with the help of the tight-binding approximation 
valid for the model of weakly coupled oscillators. This 
approach leads to the effective discrete NLS equation, 
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FIG. 8. Development of the instability-induced pe- 
riod-doubling modulations. Initial profile corresponds to a 
slightly perturbed staggered solution with Io — 29.87. The 
parameters are the same as in Fig. (H 
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C. Bright Spatial Solitons 



v(x) and w(x): 



1. Odd and even localized modes 

Stationary localized modes in the form of discrete 
bright solitons can exist with the propagation constant 
inside the band gaps, when \rf\ > 2. Additionally, it has 
been found that such solutions exist only for r/x < 0. It 
follows from Eq. (|l2|), that £ > in the IR gap and the 
first BR gaps (see Fig. ||), so that the type of the nonlin- 
ear response is fixed by the medium characteristics, since 
X = sign7. Therefore, the self- focusing nonlinearity can 
support bright solitons in the IR gap (where rj < —2), 
i.e. in the conventional wave-guiding regime. In the case 
of the self-defocusing response, bright solitons can exist 
in the first BR gap, owing to the fact that the sign of 
the effective diffraction is inverted (rj > 2). In the latter 
case, the mode localization occurs in the so-called anti- 
waveguiding regime. 

Let us now consider the properties of two basic types 
of the localized modes: odd, centered at a nonlinear thin- 
film waveguide, and even, centered between the neigh- 
boring waveguides, so that U\ n \ — x s U-\ n \-s, where 
s = 0,1, respectively. For discrete lattices, such solu- 
tions have already been studied in the literature (see, 
e.g., Ref. |^|), and it has been found that the mode pro- 
file is "unstaggered" (i.e. U n > 0) if rj < —2. On the 
other hand, Eq. (|ll| ) possesses a symmetry transforma- 
tion U n — * (—l) n U n , rj — > —rj, and x ~ * — X> which means 
that the solutions become "staggered" at 77 > 2. 
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FIG. 9. Top: power vs. propagation constant for odd 
(black) and even (gray) localized modes in a self-focusing 
(7 = +1) regime: solid — stable, dashed — unstable, and 
dotted — oscillatory unstable. Bottom: the profiles of the 
localized modes corresponding to the marked points (a,b) in 
the top plot. The lattice parameters are the same as in Fig. ^ 

To study linear stability of the localized modes, we con- 
sider the evolution of small-amplitude perturbations in 
the form (||) and, using the original model of Eqs. (||) 
and fli~C|), obtain the linear eigenmode problem for small 



+ [(a + 2~fu n )v + ju n w] S(x — hn) = 0, 

n 

+ Y~] [{a + 2ju n )w + ju n v] S(x — hn) = 0, 

n 

whe re T is related to the instability growth rate (see 
Sec. II B above). 

Our analysis reveals that even modes are always un- 
stable with respect to a translational shift along the x 
axis. On the other hand, odd modes are always stable 
in the self-focusing regime [see Fig. [)), but can exhibit 
oscillatory instabilities in the self-defocusing case when 
their power exceeds a certain critical value [see Fig. [Tcj ] . 




FIG. 10. Top: power vs. propagation constant in the 
self-defocusing (7 = —1) regime. Notations are the same as 
in Fig. B. 



2. Soliton bound states 



"twisted" modes 



Due to a periodic modulation of the medium refractive 
index, solitons can form bound states pTJI . In particular, 
the so-called "twisted" localized mode plfl i s a combina- 
tion of two out-of-phase bright solitons plfl . Such solu- 
tions do not have their continuous counterparts, and they 
can only exist when the discreteness effects are strong, i.e. 
for [77! > ?7cr- The properties of the twisted modes depend 
on the separation between the modes forming a bound 
state. We consider the cases of two lowest-order solution 
families (i) of "even" type with zero nodes (m = 0) in- 
between the peaks, and (ii) of "odd" type with one node 
(to = 1) in the middle, with Uq = 0. Then, the symmetry 
property is U\ n \ +m = — x m+1 f7_i„|_i, and we find that 
rj cr {m = 0) ~ 3.32 and rj cr (m = 1) ~ 2.95. 

Stability properties of the twisted modes in the IR 
gap, under the self-focusing conditions, can be similar 



7 



to those earlier identified in the framework of a discrete 
NLS model |^2[ . In the example shown in Fig. [n], the 
modes are stable at larger values of the propagation con- 
stant, and they become oscillatory unstable closer to the 
boundary of the existence region. Quite importantly, we 
see that the stability region is much wider in the case of 
the odd twisted modes due to a larger separation between 
the individual solitons of the bound state. 




x+h/2 x 

FIG. 11. Top: power vs. propagation constant for 
odd (black) and even (gray) twisted localized waves in a 
self- focusing (7 = +1) regime. Notations are the same as 
in Fig. |. 

The characteristics of the twisted modes in the BR gap 
can differ substantially from the previous case. First, 
the value of 77 is limited from above (2 < 7/ < 77 max ) 
and, therefore, some families of the twisted modes with 
to < TO cr may not exist. For example, for the medium 
parameters corresponding to Fig. ||(a), we have r] cr (m — 
1) < (f?max — 3.18) < ?7 cr (TO = 0), so that m CI = 1. 
Under these conditions, the even modes with to = can- 
not exist in the BR regime. On the other hand, the odd 
modes with to = 1 can exist, and they are stable in a 
wide parameter region, see Fig. [12. 
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FIG. 12. Top: power vs. propagation constant for the odd 
twisted localized waves in a self-defocusing (7 = —1) regime. 
Notations are the same as in Fig. [Tl| 



D. Dark Spatial Solitons 

Similar to the continuum NLS equation with self- 
defocusing nonlinearity |^3| or a discrete NLS model, our 
model supports localized modes on a Bloch-wave back- 
ground, also called dark solitons. Such stationary local- 
ized modes appear at the band-gap edge where 77 = 2% 
(hereafter we consider only the case of a background cor- 
responding to B loch-w ave solutions with K = 0, 7r in- 
troduced in Sec. [II B ). As the background intensity in- 
creases, the band-gap edge shifts to match the dark soli- 
ton propagation constant (3, and it follows from Eqs. jl^ ) 
and ( |I3| ) that r}(0; a) = 2\, where a = (a + 7«^ )- Since 
these solutions correspond to the first transmission band, 
we have £(/3) > (see Fig. ||), and it follows from Eq. ( p^ ) 
that the propagation constant (3 increases at higher in- 
tensities in a self- focusing medium (7 > 0) , and decreases 
in a self-defocusing medium (see Figs. [l3| and |l4|, top). 
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FIG. 13. Top: complementary power vs. propagation con- 
stant for odd (black) and even (gray) dark localized solitons 
in a self-focusing (7 = +1) regime. Notations are the same as 
in Fig. |. 

Similar to the case of bright solitons discussed above, 
two basic types of dark spatial solitons can be identi- 
fied as well, namely, odd localized modes centered at a 
nonlinear thin-film waveguide, and even localized modes 
centered between the neighboring thin-film waveguides. 
Thus, all modes satisfy the symmetry condition, U\ n \ +S — 
— (— X) s+1 U-\ n \-i, where s = 0,1, respectively. The 
background wave is unstaggered if x = — 1; an d it is 
staggered for \ = +1; the corresponding solutions can 
be constructed with the help of a symmetry transforma- 
tion U n — > (— l) n U n . However, the stability properties 
of these two types of localized states can be quite dif- 
ferent. Indeed, it has been demonstrated in Sec. Ill B 



that the staggered background in a self-focusing medium 
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(x = + 1) can become unstable. On the contrary, the 
unstaggered background is always stable if % = — 1. 
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FIG. 14. Top: complementary power vs. propagation con- 
stant in the self-defocusing (7 = —1) regime. Notations are 
the same as in Fig. [la. 



applicable for shallow modulations and a narrow gap. We 
believe our results may be useful for other fields, such as 
the nonlinear dynamics of the Bose-Einstein condensates 



ACKNOWLEDGMENTS 



The authors are indebted to Dr. O. Bang, Prof. P. L. 
Christiansen, and Prof. C. M. Soukoulis for collaboration 
at the initial stage of this project, and to Prof. A. A. 
Maradudin and Prof. Y. Silberberg for useful and en- 
couraging discussions. This work has been supported by 
the Performance and Planning Fund of the Institute of 
Advanced Studies, the Australian National University, 
and by the Australian Photonics Cooperative Research 
Center. 



Two types of the dark spatial solitons existing in our 
model are shown for both staggered and unstaggered 
background waves in Figs. |l3| and O, respectively. We 
characterize the family of dark solitons by the compli- 
mentary power defined as 

/+nh 
(\u(x + 2nh)\ 2 - \u(x)\ 2 ) dx, 
-nh 

where n is integer. The localized solutions shown in 
Fig. |l4| are similar to those found earlier in Ref. in 
the context of the superflow structures in a periodic po- 
tential. 



IV. CONCLUSION 

In the framework of a simplified model of a peri- 
odic layered medium, we have analyzed nonlinear guided 
waves and (bright, dark, and "twisted") discrete spatial 
solitons of two types, i.e. the nonlinear guided waves 
localized due to the total internal reflection and Bragg- 
like localized gap modes. We have considered two spe- 
cific geometries of the nonlinear periodic structures which 
may find their application in experimental realizations, 
namely, a single nonlinear layer embedded into a periodic 
linear medium, and a periodic array of identical nonlinear 
layers, the so-called one-dimensional nonlinear photonic 
crystal. We have analyzed the existence and stability of 
the nonlinear localized modes in both the models, and 
have described also modulational instability of homoge- 
neous states in a periodic structure of the nonlinear lay- 
ers. In particular, we have discussed both similarity and 
difference with the models described by the discrete NLS 
equation, derived in the frequently used tight-binding ap- 
proximation, and the results of the coupled-mode theory 
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